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I discuss the use of path integrals to study strong-interaction physics from first 
principles. The underlying theory is cast into path integrals which are evaluated 
numerically using Monte Carlo methods on a space-time lattice. Examples are 
given on progress related to nuclear physics. 
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1. Introduction 

Quantum chromodynamics (or QCD) is the underlying theory of the strong 
interaction (one of the four fundamental interactions in nature besides grav- 
ity, electromagnetism and the weak interaction). It is the force that binds 
quarks and gluons as the nucleus in the heart of the atom. Since all matter 
in the universe is built that way, unraveling the structure of matter at its 
deepest level as governed by QCD is key to our understanding of the phys- 
ical world, and presents one of the most challenging tasks in contemporary 
nuclear and particle physics. While it is easy to write down the basic equa- 
tion of QCD, it is very difficult to obtain quantitative solutions in view 
of the complex quark-gluon dynamics. At present, the only known way to 
solve QCD directly is by numerically evaluating the path integrals in the 
theory on a discrete space-time lattice using supercomputers. For textbooks 
on the subject, see Refs.^"^ 



2. Formalism 

All physics can be computed by path integrals in QCD. Take the calcula- 
tion of the proton mass for example. It requires the fully-interacting quark 
propagator defined by the Euclidean-space path integrals over gluon field 
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and quarks fields ip and tp, 

The action is the sum of a ghion part and a quark part 

Sqcd = SG + Sg = ^J dx^TrF^, F"^ + J dx^i^M^P, (2) 

where F^i, = d^A,y — +g[Af^, A^,] is the gluon field strength tensor and 
M = j'^Df^ + TUq is the Dirac operator for quarks and 7^ tlic^ 4x4 gamma 
matrices. The interaction between the two is provided by the co variant 
derivative = + gA^. The basic parameters of the theory are the 
coupling constant g and quark mass niq. The quark part of the integration 
can be done analytically using Grassmann variable integration, leading to 
a path integral over only gluon fields 

1 f DAadct(M)e-'^'=^M-^ 

which is evaluated numerically by Monte-Carlo methods. In the quenched 
approximation, the quark determinant det(M) is ignored (equivalent to 
suppressing vacuum polarization), making the numerical calculation signif- 
icantly faster. 

To compute the proton mass, one considers the time-ordered, two-point 
correlation function in the QCD vacuum, projected to zero momentum: 

G{t) = J2{0\T{rj{x)m}\0). (4) 

Here rj is called the interpolating field which is built from quark fields with 
the quantum numbers of the proton (spin-1/2, isospin-1/2, quark content 
uud, charge -|-e) 

r]{x) = e"^" [u'''^{x)C^5d\x)] u%x) (5) 

where C is the charge conjugation operator and the superscript T means 
transpose. Sum over color indeces a,b,c is implied and the e"'^" ensures the 
proton is color-singlet. 

The calculation of G{t) at the quark level proceeds by contracting out 
all the quark pairs, resulting in 

Git) = Y^e'^bc^'^'"'^' {sfj,CSf^Cj,Sf +SfTr{GS^'^Zci5St''j5)} , 

X 

(6) 
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On the hadronic level, the correlation function is saturated by the com- 
plete spectrum of intermediate states (after a parity projection) with the 
proton as the ground state: 



where rrii are the masses and are the 'amplitudes' which are a measure 
of the ability of the interpolating field to excite or annhilate the states 
from the QCD vacuum. At large time, the ground state proton, e""*!*, 
dominates G{t), with the excited states exponentially suppressed. Other 
physics quantities are computed in a similar way. 

3. Some examples 

It is impossible to give; a full account of the achievements made in lattice 
QCD in the space given here. A good place to gauge the progress in the 
field (which is not limited to QCD) is the annual gathering'* by active 
practitioners. Here I select a few examples relevant to nuclear physics. 

Example 1. Fig. 1 shows that the computed light hadron mass spec- 
trum which comes within 10% of the observed spectrum. The remaining 
discrepancy is attributed to the quenched approximation. Intense eff'orts 
are under way to extend the success to the excited sectors of the mass 
spectrum.^ 

Example 2. Fig. 2 shows the high-precision lattice QCD calculations 
for selected physics quantities which are within 3% of the observed values. 
The effects of the quenched approximation are clearly demonstrated. 

Exemiple 3. Fig. 3 shows a calculation of the magnetic moment of the 
proton ans neutron which find good agreement between lattice QCD and 
experiment. Results from a different method can be found in Ref.,^ as well 
as precise determinations of the strangeness magnetic (Gm)^*' and electric 
form factors (G|;)*^ of the proton. 

Example 4. Fig. 4 shows the good agreement between lattice QCD 
and experiment for pion-pion scattering length,*^ using Liischer's method. ""^^ 
Other hadron-hadron scattering channels (7r+if +, nA, nS~, etc.) have been 
studied by the same group. 

Example 5. Fig. 5 shows a direct lattice QCD calculation of the 
nuclcon-nuclcon potential*"' whose features arc consistent with the known 
phenomenological features of the nuclear force. 
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Fig. 1. The computed light hadron mass spectrum in the quenched approximation. See 
Ref.^ for details. 
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Fig. 2. Lattice QCD results compared to experiment as ratios for selected physics quan- 
tities (the dashed line is perfect agreement), in the quenched approximation of QCD (left 
panel) and full QCD (right panel). See Ref.^ for details. 
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Fig. 3. Magnetic moment of the proton and neutron as a function of the pion mass 
squared. The hne is a chiral extrapolation. See Ref.* for details. 
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Fig. 4. Lattice QCD calculation of the tt — tt scattering length in the isospin 2 channel. 
The dashed line is a chiral fit. See Ref.^^ for details. 



4. Conclusion 

The path-integral formulation of QCD, coupled with large-scale numerical 
simulations, has played a crucial role in the progress of nuclear and particle 
physics. Its continued prominence as the only way to solve QCD with con- 
trolled errors is expected to play out with better results in the foreseeable 
future. This work is supported in part by U.S. Department of Energy under 
grant DE-FG02-95ER40907. 
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Fig. 5. The lattice QCD result of the central (effective central) part of the NN potential 
in the ^Sq (^Si) channels. The inset shows its enlargement. The solid lines correspond 
to the one-pion exchange potential (Yukawa potential). See Ref.^* for details. 
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